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Abstract 

The multifractal properties of electronic eigenstates at the metal-insulator 
transition of a two-dimensional disordered tight-binding model with spin-orbit 
interaction are investigated numerically. The correlation dimensions of the 
spectral measure D 2 and of the fractal eigenstate D 2 are calculated and shown 
to be related by D 2 = 2D 2 . The exponent 77 = 0.35 ± 0.05 describing the 
energy correlations of the critical eigenstates is found to satisfy the relation 
rj = 2-D 2 . 

PACS numbers: 71.30.+h, 64.60.Ak, 72.15.Rn 
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Electronic states in disordered systems are known to be either localized or extended. 
At T = the so called mobility edge separates insulating (localized) from current carrying 
(extended) states. Recently, the electronic properties directly at the critical point have 
received increasing attention [|l]-|J. Two-dimensional (2d) systems are very well suited for 
numerical investigations of the electronic eigenstates at a critical point. The only systems, 
however, that exhibit a complete Anderson transition in 2d are models with symplectic 
symmetry. This is in contrast to the orthogonal and the unitary case where all states are 
localized (weak localization) and the quantum Hall (QHE) systems where only localized 
states with a diverging localization length at some singular energies can be observed. 

In a recent paper, Chalker and coworkers |7j reported on numerical investigations of eigen- 
state fluctuations and correlations near the mobility edge of a two-dimensional tight-binding 
model with spin-orbit coupling. It was found that the probability amplitude distribution 
exhibits multifractal behavior which can be characterized by a set of generalized dimensions 
D q . The fractal spatial structure of the wavefunctions also shows up in the eigenfunction cor- 
relations between states close in energy. A similar behavior was found also in the quantum 
Hall model 0,§-(TC], i.e. a two-dimensional disordered system of non-interacting electrons 



where a strong magnetic field causes the localization length to diverge at the centres of 
the Landau bands with a universal exponent For these critical states, characteristic 



quantities like the generalized correlation dimension of the wavefunction D 2 = 1.62 |8]|J, the 
corresponding generalized dimension of the spectral measure D 2 = 0.81 @, the exponent 
rj = 0.38 [|l]|| governing the energy eigenfunction correlations, and «o = 2.3 P,|10| which 



for — 1 < q < 1 determines the so called /(o;)-distribution are known, partly with sufficient 
precision. This is, however, not the case for the spin-orbit system mentioned above, mainly 
because of the small system sizes considered so far (M < 18 a) 0. The authors of Ref. 
used a model proposed by Evangelou and Ziman |T2[ which represents a two-dimensional 
disordered electronic system with spin-orbit interaction. 

In this letter, a different lattice model with symplectic symmetry is investigated for 
system sizes up to 150 a x 150 a, where a is the lattice constant. It is found that D 2 = 
1.66 ± 0.05, D 2 = 0.83 ± 0.03, r) = 0.35 ± 0.05, and a = 2.19 ± 0.03. Thus, within 
the numerical uncertainties, the proposed relations rj = 2 — D 2 |T3],3, D 2 = 2D 2 ||, and 
Ag yp 7r = (ctQ — 2) _1 [ 14 1 hold also for the two-dimensional system with spin-orbit interaction. 



The model used to calculate the eigenstates has been put forward by Ando [p~5| , p~6|| to 
simulate two-dimensional systems in n-channel inversion layers on surfaces of III- VI semi- 
conductors. The Hamiltonian describing this situation is 

W = £ m C L,aCm,a' + £ V ^ ^ U > a ') C ln,* C n,v', W 
m,cr mn,cnj' 

with disorder potentials e m , creation $ ma and annihilation c ma operators of a particle 
at site m and spin state a, respectively. The transfer matrix elements V(m, cr; n, a') = 
a a' V(n, —a'; in, —a) which are restricted to nearest neighbours only depend on whether the 
transfer from site m to one of the nearest neighbours n goes along the x- or the y-direction. 
The strength of the spin-orbit interaction is determined by the parameter S = V 2 /V, with 
V = (V? + V 2 2 ) 1/2 taken to be the unit of energy, where V\ and V 2 are the matrix elements 
for transitions with and without spin-flip, respectively. V(m, a; na') is then given by 
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y-= ( r \, 1, v v =(v\ JM, k = +i) = i : i, k = -i> = ( : i. ( 2) 
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The localization properties of symplectic models have previously been analyzed numerically 



,I3-P0[ from which a metal-insulator transition can be inferred. The most recent calcu- 
lations of Fastenrath ||20|j2l|| report a critical exponent v = 2.75 for the localization length 
at the band centre, E/V = 0, together with a critical disorder W c = 5.74V for a constant 
probability distribution of the on-site disorder potentials {e m } and a spin-orbit strength 
S = 0.5. We also take these parameters in what follows so that our new results supplement 
the already published data. The eigenvalues and eigenstates were calculated numerically by 
means of a Lanczos-algorithm for systems of size up to L = 150 a with periodic boundary 
conditions applied in both directions. 



The structure of the eigenstates is analyzed in terms of the /(a)-distribution p2| , |23|| which 
completely characterizes the spatial scaling behavior of the g-moments of the wavefunction. 
In Fig. 1 the /(a)-distribution for a particular eigenstate from the critical region near the 
center of the tight binding band (E = V) and a disorder strength W = 5.74 V is shown for 
several g-values. For comparison the corresponding parabolic approximation [0,0 i s a l so 
shown. The average over 190 eigenstates taken from the energy interval [—0.15,0.0] gives a 
value of a = 2.19 ± 0.03. This number represents the most probable value of the scaling 
exponents a q = d/dq (q — l)D g . 

Using the box-probability method the correlation dimension D2 of the multifractal eigen- 
states is obtained from the scaling of the second moment (q=2) of the averaged box- 
probability P(q, A) = Ei(ErefMA) l^( r )l 2 ) 9 ~ \ {q ~ 1)D \ where f^(A) is the size Ixl of the i-th 
box. A power law relation is observed for / = XL in the range 2a < I < L/2. The average 
over all eigenstates from the energy interval gives a correlation dimension (-D 2 ) e = 1.66±0.05 
which is close to the value of 1.63 reported in for smaller system sizes. 

The correlation dimension of the spectral measure, D2, which is related to the temporal 
decay |25||| of the maximum of a wavepacket (probability of return) built from critical 



eigenstates was calculated for the same energy interval. This exponent is obtained from the 
scaling relation of the local density of states 



7 (g,.) = liml/L 2 £E( £ \a E \^ ^ e (*-D*,, (3) 

with dE = ^e{ v )/^jE' |'0-E'( r )| 2 ) 1/ ' 2 - The scaling behavior is shown in Fig. 2 for q = 2. 
Here, an exponent D 2 = 0.83 ± 0.03 is found which satisfies the relation D 2 = 2D 2 proposed 
previously for the 2d QHE system |J. These results indicate that the diffusion at the 
mobility edge will be non-Gaussian in the long-time limit. Due to the multifractal nature of 
the spatial amplitude fluctuations of the critical wavefunctions and the behavior of the fractal 
spectral measure, the diffusion coefficient will not be a constant. A non-trivial frequency 
and wavevector dependence governed by an exponent 77 has been observed in the QHE model 
Hill] and a similar behavior was also seen for the symplectic case 0. 
In Fig. 3 the correlations of the eigenstates close in energy 

Z(E, E>) = ]T |tMr)| 2 hMr)| 2 ~ \E - E'\^ d (4) 
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averaged over small energy intervals are shown for the critical eigenstates. A power law 
relation is observed with 77 = 0.35 ± 0.05. From the above results it is seen that the relation 
7] = d — D 2 [p~3|,|9|j also holds in the 2d symplectic case. 



We have to mention that our value for 77 is compatible with the result for a different 
symplectic model obtained from the calculation of the two-particle spectral function at 
the mobility edge pj. It is also in accordance with an earlier estimate by Evangelou J26 



although the corresponding scaling exponent of the localization length, v = 1.6, differs 



considerably from v = 2.75 obtained by Fastenrath [^Ol^TJ for the Ando Model. However, 
to determine rj from the calculated finite size scaling variable at the critical point, A c , 
Evangelou used the relation A c = 1/ (tti]/2) which is only within the parabolic approximation 
(D q = 2 — q(a — 2) |24l , |l4ll , valid for \q\ < 1) equivalent to a relation proposed recently 
by Janfien A* yp = l/(ir(ao — d)), relating the typical finite size scaling variable at the 
critical point, A* yp , to the multifractal behavior of the eigenstates. 

We note that our results for a , D 2 , D 2 , and 77 are very close to those obtained for the 



2d QHE model |T0|,|3|. Although we could not observe any size dependence, the presently 



achieved system sizes do, however, not allow to exclude the possibility that the values 
actually coincide in the thermodynamic limit, L —>■ 00. A similar correspondence has recently 
been asserted for the energy level statistics || of the critical 3d Anderson model with and 
without a magnetic field. In addition, the critical exponent of the localization length was 



reported to be identical in both cases |27||28| . If these observations for 3d models were correct 
and if also the values obtained from the multifractal analysis in 2d for the quantum Hall 
systems were indeed the same as those reported above for the symplectic disordered systems, 
then the behavior at the critical point would primarily be determined by the euclidean 
dimension and not by the symmetry class of the Hamiltonian. 

In conclusion, the multifractal properties of the electronic eigenstates at the metal- 
insulator transition of two-dimensional disordered systems with symplectic symmetry have 
been investigated. The /(a)-distribution, the correlation dimensions of the spectral mea- 
sure and the critical eigenstates, and the energy correlations of the wavefunctions were 
calculated. The obtained values for a = 2.19 ± 0.03, D 2 = 0.83 ± 0.03, D 2 = 1.66 ± 0.05, 
and 7] = 0.35 ± 0.05 appear to be independent of system size and satisfy general relations 
that have been proposed previously for the critical states at the metal-insulator transition 
in other disordered systems. 
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FIGURES 



FIG. 1. The /(a ((^-distribution function of critical eigenstates calculated for a system of size 
L/a = 150 and q = ±3, ±2, ±1.5, ±1, ±0.8, ±0.5, ±0.3, 0. The full curve is a fit using the 
parabolic approximation with = 2.19. 

FIG. 2. The scaling of the spectral measure at the mobility edge for a system of size L/a = 100 
from which a correlation exponent D2 = 0.83 ± 0.03 is obtained. 

FIG. 3. The energy correlation Z(E, E') of the critical eigenstates as a function of energy 
separation showing a power law relation ~ \E — E'\^ v ^ 2 with an exponent 77 = 0.35 ± 0.05. 



6 



